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Abstract. Binary galaxies are modeled as point-masses obeying the non-Newtonian 
MOND and Mannheim-Kazanas (MKG) theories of gravity. Random samples of such 
systems are generated by means of Monte Carlo simulations of binary orbits. Model pairs 



I/") have total masses and mass ratios similar to pairs in the cataloged sample used in the 

o 

analysis. 



ON 



General features of synthetic samples are derived from a comparison with observed 



O-ndata of galaxy pairs in R x AV/(L\ +L2) 1 / 2 space. Both MOND and Mannheim-Kazanas 
2 binaries either on circular or low-eccentricity orbits cannot be the source of observations 
^ because they require extremely high M/L values (ps 45 solar units). Both MOND and 
MKG binaries on high-eccentricity orbits and reasonable M/L values (5 solar units) 
produce envelopes of R x AV/{Li + L'z) 1 ^ 2 consistent with the observations, but the 
distribution of separations is inconsistent with the observed data, unless strong selection 
effects are at work. 

A definite answer to the issue whether one or another model is suitable to explain real 
binary galaxy dynamics will be only possible when a large sample containing a significant 
fraction of wide pairs, determined with velocity -blind selection procedures, is investigated 
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contamination by non-physical pairs. 
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1. Introduction 

The dynamics pertaining to Newtonian gravity are unable of avoiding mass discrepancies 
in astronomical systems unless some kind of, yet unknown, dark particles are invoked 
(Faber and Gallagher 1979, van Albada and Sancisi 1986, Trimble 1987, Sanders 1990, 
Ashman 1992). For example, the MACHO and EROS searches for dark matter in the 
form of compact sub-luminous objects in our Galaxy halo (Alcock et al. 1993, Bennett 
et al. 1995, Aubourg et al. 1995, and references therein) have not detected sufficient 
microlensing events to explain the kinematics of the Milky Way with a spheroidal halo 
populated by brown dwarfs. Furthermore, direct optical searches of the Hubble Space 
Telescope have not also been able to detect the predicted amount of matter in the form of 
faint stars in the halo of the Milky Way (Bahcall 1994). Another much explored approach 
to the problem has been the suggestion that Newton laws are not adequate and that the 
solution might be found in alternatives to classical gravity and dynamical laws. Many of 
such alternatives are empirically motivated (e.g., Milgrom 1983a, 1983b, 1983c, Sanders 
1984, Kuhn and Kruglyak 1987) by the observations of flat rotation curves but that is 
not always the case. Mannheim and Kazanas (1989, 1991) have derived an alternative 
theory of gravity whose main motivation was the search for the correct covariant general 
theory of gravity. Specifically, they found the exterior solution to conformal Weyl gravity 
associated with a static, spherically symmetric gravitational source. 

The most successful of the empirical alternatives to Newton's laws has been Milgrom's 
MOND (Modified Newtonian Dynamics) which has been subjected to successive experi- 
mental tests (e.g., Gerhard 1994, but see Milgrom 1995). The theory has been strongly 
and successfully defended by Milgrom in many papers in recent years (see also Bege- 
man, Broeils and Sanders 1991, Sanders 1994). Nevertheless, a most serious setback is 
that there is no covariant theory of gravity that, in the weak field limit, reduces itself to 
MOND. Some attempts have been made otherwise in this direction (Bekenstein and Mil- 
grom 1984, Sanders 1986) but yet MOND remains as a phenomenological ad hoc theory, 
being this, up to now, its main source of criticism. Other consequences, dynamical and 
cosmological, of MOND are discussed by Felten (1984). 



(hereafter MKG) to model the circular velocity profiles of four spiral galaxies with fairly 
reasonable fits, without the requirement of dark matter. He claims that Einstein-Newton 
gravity proves to be inadequate in comparison with fourth order conformal Weyl gravity, 
represented by the MKG solution. In spite of that, and of having on its foundation a 
general theory of gravity with the status of a full covariant one, MKG has been ques- 
tioned on its ability of explaining exceptionally extended flat rotation curves (Sanders 
and Begeman 1994). Moreover, even conformal Weyl gravity has been criticized on cos- 
mological grounds. Elizondo and Yepes (1994) presented exact solutions to the conformal 
Weyl gravity cosmological equations that failed to yield primordial nucleosynthesis abun- 
dances consistent with present observational constraints. They conclude that conformal 
cosmological models are very unlikely to give a realistic description of the Universe. 

This is intended to be the first one of a series of papers devoted to binary galaxy 
dynamics, under the premises of alternative theories of gravity. Here, it is investigated 
what binary galaxy data might suggest one about MOND and MKG. Simple point- 
mass models are used to represent bound pairs of galaxies under MOND and MKG 
prescriptions for the gravitational mutual galaxy interaction. Monte Carlo simulations 
of synthetic samples are performed by means of numerical calculations of pair orbits. 
They are subsequently projected on the sky and the relevant observed quantities derived. 
A qualitative comparison between a sample of real pairs and simulated samples is then 
made. The second paper of the series will present a rigorous statistical analysis of observed 
samples of binaries in the light of MOND, and the third one will focus on MKG. Biases 
introduced by selection effects on sample determinations and "contamination" by non- 
physical pairs must be considered separately for each theory of gravity. It is obvious, for 
example, that the definition of what a non-physical pair is depends upon the theory that 
describes the interactions and the dynamics of galaxies in pairs. 

Section 2 presents the data which come as a list of disk binary galaxies extracted 
from the Catalogue of Multiple Galaxies. A fiducial Keplerian model, determined in the 
Appendix, is qualitatively fitted to the data. Such a model is meant to be confronted with 
MOND and MKG simulated samples. In section 3, the MOND gravitational potential for 
point-mass binaries and technical details of the simulations are described. The same is 
done in section 4 for MKG. Section 5 is devoted to Monte Carlo simulations of MOND 
and MKG synthetic samples. In section 6, results of the qualitative comparison between 
observed and synthetic binaries are discussed and the main conclusions presented. 



Much effort has been dedicated in the last years to the compilation of binary galaxy 
lists (Karachentsev 1972, 1987, Turner 1976, Peterson 1979, Schweizer 1987, Soares 1989, 
Chengalur, Salpeter, Terzian 1993, 1994, Soares et al. 1995, etc). For the sake of unifor- 
mity, only one list of pairs is considered, namely, that extracted from the Catalogue of 
Multiple Galaxies (CMG, van Moorsel 1982, 1987, Oosterloo 1988, 1993, Soares 1989). 
The selection of pairs was done by applying a surface density enhancement procedure 
to the Uppsala General Catalogue of Galaxies (Nilson 1973), a method devised by T.S. 
van Albada. A description of the method can be found in van Moorsel (1982), Soares 
(1989, chapter 2, by van Albada and Soares) and Soares et al. (1995). The method does 
not require any redshift information, and this is a feature that makes it particularly 
suitable for investigations concerning alternative physics. The concept of a bound pair is 
dependent on what gravity law rules the dynamics of the pair, and cannot be relied upon 
simple upper limits on line-of-sight velocity differences, which is a typical characteristic 
of velocity-dependent selection procedures. Still, on uniformity, from the whole list of 
358 binaries in the CMG a homogeneous subsample of 230 pairs is considered (see Soares 
1989), which is formed by rotationally supported galaxies, ranging from lenticulars to late 
spirals. The main reason for avoiding pairs with elliptical galaxies is that both MOND 
and MKG have gathered their status of likely alternative theories after being applied on 
modeling circular velocity radial profiles of spiral galaxies. The whole sample is available 
on machine readable form upon request (dsoaresfif isica.ufmg.br). The line-of-sight 
heliocentric velocity differences are normalized by the square root of the total blue lu- 
minosities of the pairs, in units of 10 11 Lq. Such a normalization allows one to scale the 
envelopes given by eqs. All and A14 (Figure 6 — Appendix) by a proper M/L value. 

It has been shown in many recent binary galaxy investigations, based on conventional 
physics, that there might be an extended dark halo either around the individual galaxies 
or embedding the whole pair (White et al. 1983, Schweizer 1987, Soares 1989, Charlton 
and Salpeter 1991, Bartlett and Charlton 1995). Although the presence of optical (i.e., 
non-physical) pairs in a sample could bias the analysis against point-mass models (as 
may have happened with the study by White et al. 1983, as pointed out by Pichio and 
Tanzella-Nitti 1985, and Pacheco and Junqueira 1988) there is compelling evidence that 
Keplerian models are far from being a good description of binary dynamics. Nevertheless, 
Keplerian approximations serve as reliable upper limits to AV/ L 1 ! 2 for a sample of binary 
galaxies, provided one not be interested in the detailed distribution of pairs in that space. 



ies in pairs as a function of projected separation in the plane of sky are derived. Such 
envelopes are here confronted with real data and set definite limits upon the overall dis- 
tribution of the observations. It must be pointed out that the Keplerian model derived 
in the Appendix serves only as reference. 
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Fig. 1. Keplerian envelopes in the R — AV/L 1 ^ 2 plane (curves). Filled circles represent CMG 
pairs. There are three model free parameters: the mass-to-light ratio (M/L), the minimum 
pericentric separation (r m ) and the orbital eccentricity (e). The solid curve, in the four panels, 
is the fiducial envelope mentioned in the text, with M/L = 30, r m = 20 kpc and e = 0.4, chosen 
as a reference frame for CMG pairs [in = (L\ + L2) /^-O 11 Lq\. 

Figure 1 shows how to fill in the binary phase space with varying the three free parameters, 
M/L, r m (minimum pericentric separation) and e (orbital eccentricity), allowed by the 
Keplerian envelope models derived in the Appendix. The ordinate AV / L 1 ! 2 is obtained 
by multiplying AV/M 1 ! 2 by the scaling factor (M / L) 1 ! 2 . The bottom right panel shows 
the CMG pair sample. 

The solid curve in Fig. 1 was chosen as a fiducial envelope for the subsequent analysis. 
Its parameters are M/L = 30, r m = 20 kpc and e = 0.4. From a Newtonian point of 
view, pairs above it are most likely unbound pairs, the so-called optical pairs, because 



and MKG as it is shown in the Monte Carlo simulations below. 




Fig. 2. The envelope of maximum AV/L 1 ^ 2 is shown as a histogram, determined after exclusion 
of outliers. The number of pairs in each bin and the clipping weight are given in the right bottom 
corner of the panel. The smooth curve is the fiducial Keplerian envelope. 

To verify whether such a Keplerian model is a fair reference envelope for the observations, 
a ka clipping procedure (see Yahil and Vidal 1977) was applied to the sample to exclude 
outliers, which would bias the determination of the envelope of maximum velocities. 
Figure 2 shows the resulting envelope confronted with the fiducial Keplerian envelope 
chosen above. Each bin in the histogram of Fig. 2 has 50 pairs. For every bin, the mean 
and standard deviation of AV/L 1 / 2 are calculated, and pairs with AV/ L 1 ! 2 over 2.5a 
are excluded from the initial bin population. A new mean and a are calculated and 
the process is repeated. The iteration ends when every pair falls within 2.5u. Then, the 
highest AV / L 1 ! 2 is computed. It is apparent from Fig. 2 that the chosen parameters for 
the Keplerian envelope are consistent with the numerical envelope derived from the data. 



The most successful phenomenological theory, alternative to dark matter Newtonian ap- 
proach to the mass discrepancy problem, has been put forward by Milgrom, originally in 
three successive papers (1983a, 1983b, 1983c). It consists of a modification in Newtonian 
dynamics in the regime of low accelerations. In the framework of Milgrom's Modified 
Newtonian Dynamics (MOND), galaxies have no dark mass. Its simplified formulation 
has two basic assumptions: (i) Newtonian dynamics break down when the accelerations 
involved are small, and (ii) the acceleration a of a test particle at a distance r from a 
mass m is given by a 2 ja k, Gmr~ 2 in the limit Gmr~ 2 <C o, (or a <C £*<>)• The constant 
a plays the role of a transition acceleration (implying a transition length scale, r t , which 
is a function of the mass m) from the Newtonian to the MOND regime. In other words, 
in regions where the acceleration is smaller than the limit acceleration the gravitational 
field is roughly proportional to the inverse of the distance from the center of the galaxy, 
and as a consequence, spiral galaxy flat rotation curves can be obtained. The value of a 
was originally determined by Milgrom (1983b) to be k, 4 x l(T 8 cm/s (A Hubble constant 
H = 75 kms^Mpc -1 will be assumed throughout this paper). Later, a larger and better 
sample of flat rotation curves were analyzed and yielded the value of 1.21 x 10 _8 cm/s 2 , 
for the same H (Kent 1987, Milgrom 1988 and Begeman, Broeils and Sanders 1991); 
such a value is used below in the present MOND simulations of binary galaxies. 

The proper question here is in which way MOND works with binary galaxies. Milgrom 
(1986, 1994) gives numerical solutions of three problems in the framework of the Modified 
Newtonian Dynamics. They are: (1) the force law between two point masses, (2) rotation 
curves of various model disk galaxies, and (3) the field of a point mass in a constant 
external acceleration field. Of these, the first one is used here for the account of binary 
galaxy dynamics. 

Milgrom finds, in problem (1), that the MOND force field is very weakly dependent 
on the ratio M1/M2 = q of the two point masses. In the limit of large r the force F can 
be put in the form 



F(M 1 ,M 2 ,r) = A( q) , M lM2 (Ga fl 2 
K ' ' ; W (Mi +M 2 )!/2 r ^ > 

^(?) = |r 1 (i + ?) 1/2 [(i + ?) 3/2 -g 3/2 -i 

This is the force field that one expects from a logarithmic potential. In the range 1 < 
q < 6, characteristic of the CMG binary sample (see below), the dimensionless function 



The test-particle circular velocity (V ro t,i) of an individual galaxy, with mass equal to 
Mi, at distances larger than r t , is 



Vrotl = (GMiflo) 1 / 4 , (2) 

and the circular velocity [V r c ; rc (r)] of a binary system, with total mass equal to Mi + M2, 
is given by 

V clIC (r) = A(q) 1 l' 2 [G(M 1 + M 2 )a ] 1 / 4 . (3) 

Equation 3 is valid in the limit of low accelerations, or, for binary galaxy separations 
much larger than the corresponding r% = (GM / 'do) 1 / 2 , where M is the total binary mass 
Mi + M'2- In the Monte Carlo simulations, described below, the total mass of pairs will 
range from 2.9 x 10 11 Mq to 1.3 x 1O 12 M0. The transition radius for such systems ranges 
from 10 to about 20 kpc. Model galaxies have an assumed optical radius of about 10 kpc, 
thus, the binary potential is always MOND potential in the regime of r ^ r t , as long as 
the visible galaxies do not overlap. 

Each synthetic binary coming out of the simulation is obtained from the associated 
binary orbit. MOND orbits must be calculated numerically because the logarithmic po- 
tential does not admit an analytical solution. The starting point of the orbit calculation 
is the apocentric separation. Thus, one needs to know the apocentric separation r apo 
itself and the velocity at apocenter. For a given eccentricity and apocentric separation, 
the velocity is univocally defined. It is convenient to express V(r apo ) as a fraction of the 
circular velocity at apocenter, i.e., as a x V c i TC (r apo ). In the simple case of a Keplerian 
orbit with eccentricity e, the parameter a is given by 

a (e) = (l-e) 1 / 2 . (4) 

The relationship between a and the orbit eccentricity e, for the MOND binary potential, 
and galaxy separation r ^ r t , is 

where e is given by its general definition e = (r apo — r per )/(r apo -|-rp er ). The function a(e), 
valid only in the low acceleration regime, is represented in Figure 3. An artificial sample 
of MOND binary galaxies is generated in section 5 by means of Monte Carlo simulations 
of binary orbits. 
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Fig. 3. Apocentric velocity in units of apocentric circular velocity as a function of orbital 
eccentricity for binary galaxy orbits. Solid curve: MOND potential, in the low acceleration regime 
(eq. 5). Dashed curve: Keplerian potential shown as a reference curve. 

4. Mannheim- Kazanas gravity 

Mannheim and Kazanas, in a series of papers (1989, 1991, 1994; Mannheim 1992, 1993, 
1994a, 1994b) put forward the idea that gravity be based on the fourth order conformal 
Weyl theory rather than on Einstein's general relativity. While Einstein's gravity has a 
lot of success in a wide range of experimental facts, it faces its unequivocal inconsistency 
in the weak field limit, represented by Newtonian gravity, with the mass discrepancy 
problem in astronomical systems. Mannheim claims that one of the advantages of MKG 
conformal gravity is its ability, in the weak field limit, to solve the dark matter problem 
in spiral galaxies, explicitly posed by flat rotation curves. A conceptual point stressed by 
MKG's authors is that while alternative physics theories, such as MOND, depart from the 
phenomenological problem, which they are meant to explain (e.g., flat rotation curves), 
and then go on seeking a general covariant formulation, in order to get the desirable 
status of an acknowledged theory of gravity, Mannheim and Kazanas take the inverse 
way. 

Instead of being based upon a second order action, as Einstein's general relativity, 
MKG follows from a conformal invariant fourth order action which implies a fourth order 
Poisson equation (Mannheim 1994b). The latter was integrated to yield an exact exterior 
solution to a spherically symmetric gravitational source. Neglecting terms similar to the 
standard general relativistic corrections to the Newtonian potential and a high order 
term in r (see below), the following non-relativistic gravitational potential is derived, 
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where c is the speed of light and 7 a constant. The circular velocity of a binary system, 
with total mass M\ + M2, is then 



The linear potential component in eq. 6 is responsible for the good fits obtained for flat 
rotation curves (e.g., Mannheim 1993). The actual value of 7, for a typical galaxy, is 
?s 10 -28 cm -1 (Mannheim and Kazanas 1989), which "intriguingly is roughly the value 
of the inverse Hubble length", in the authors' words. In passing, Milgrom (1989) has also 
called the attention for "one fact that may prove of prime importance", i.e., "the near 
equality of a , [the MOND constant] as determined from galaxy dynamics, and that of 
cH , where H is the present value of Hubble constant". Indeed, the value of a adopted 
in this work (see section 3) is equal to 0.16 cH , and the value of 7 (see below) is 0.043 
H /c, so one derives a /j = 3.7c 2 . 

The reality of such relationships between a , 7 and H would imply that the value of 
H itself enters (non-linearly, see section 6) into the formalism of both MOND and MKG 
for interpreting binary orbits. 

The MKG constant 7 is not an universal one, like MOND's a . It must be determined 
for every single galaxy, contrary to what is found with MOND (Begeman, Broeils and 
Sanders 1991). For the four galaxies studied by Mannheim (1993), I/7 ranges from 1.3 to 
4.0 x 10 29 cm. In the simulations of section 5, eq. 6 was used with the average 7, which has 
the inverse value of 2.5 x 10 29 cm. But the fact that 7 can vary from galaxy to galaxy must 
be taken into account when a sample of binaries with a substantial fraction of wide pairs 
is investigated. It must be pointed out however that it is not clear by which mechanism a 
law of gravity may have a tunable "constant" for every gravitational system. This seems 
to be a major weakness of MKG unless there is some constraint (of cosmological origin 
or of other kind) to fix 7 as a universal constant. 

Orbits of point-masses under the potential given by eq. 6 are numerically calculated. 
The auxiliary parameter a, described in the previous section, is related to the orbital 
eccentricity in a more complicated way than in Keplerian and MOND cases. The re- 
lationship includes dependences on the total binary mass M\ + M2 and on the orbital 
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Fig. 4. Same as Figure 3 for the MKG potential (eq. 6). Solid curves: eq. 8 for r apo = 20, 50, 100 
and 400 kpc, and Mi +M2 = 1 x 10 11 Mq. Dashed curve: Keplerian potential shown as a reference 
curve. 

In deriving eq. 6, some terms were neglected, including a higher order term in r. 
Mannheim and Kazanas (1989) and Mannheim (1995) speculate that on the largest scales, 
the contributions to such terms, and even to the linear component, from all galaxies, 
might merge and add to a preexistent general cosmological background. It is difficult 
to anticipate the gravitational influence of such a cosmological background, until a full 
development of MKG implications is available. It must be pointed out, however, that one 
deals here with intermediate scales (ps 1 Mpc and less), which are probably only little 
affected by the general cosmological version of eq. 6. 

5. Monte Carlo simulations 

The Monte Carlo simulations shown here follow the general recipe given in Soares (1990, 
hereafter S90), where a study of binary galaxies with dark halos was made. In that 
study all binaries had fixed r apo = 200 kpc. Now, for the sake of generality, a power-law 
distribution of binary spatial separations is: 

P(r)dr oc dr, (9) 



radius r and external radius r + dr. The exponent Q, in the one-dimensional distribu- 
tion P(r), is approximately 0.2, if the two-point correlation function (which is usually 
expressed as a three-dimensional power-law distribution with power —1.8, e.g., Peebles 
1980) is considered to be valid on scales of binary galaxies (Gott and Turner 1979, van 
Moorsel 1982, White et al. 1983). As in S90, a lower limit in r of 20 kpc is set, dictated 
by the assumed size of a binary system at closest approach. The model binaries, both in 
MOND and MKG, have equal galaxies with 10 kpc of radius. An upper limit on r of ps 1 
Mpc, although rather arbitrary, is also chosen. It is, otherwise, consistent with the widest 
pairs in the CMG sample and with recently determined wide pair lists (e.g., Chengalur 
et al. 1993). In other words, it is implicitly assumed that pairs of galaxies separated by 
more than ps 1 Mpc are optical, be under the MOND or the MKG potential. The analysis 
made in the next section, however, is by no means affected by the choice of such an upper 
limit in r. 

For technical reasons (see below), eq. 9 is applied to binary apocentric separations 
instead of spatial separations. There is no major problem with this procedure. In case 
of circular orbits, the distributions of spatial and apocentric separations are, of course, 
coincident, and, in case of non-circular orbits they are not significantly different due to 
the fact that orbiting bodies spend most of their orbital time near apocenter. 

The simulated pairs need to have a distribution of total luminosities and luminosity 
ratios similar to the observed sample. The total masses of the simulated pairs are derived 
assuming M/L = 5, as found for MOND by Begeman, Broeils and Sanders (1991), and 
characteristic of stellar populations. Rather than adopting the distribution of observed 
(LijL'z) for the simulated pairs, which requires computing 230 orbits, the simulated 
(LijL'z) are obtained from the observed ones, averaged in a two-dimensional grid of 
(q,Li), where q = Li/L't- Table 1 shows this grid. 209 pairs (91% of the total) have q 
less than 6. Mass ratios of the order of 10 or more are typical of satellite systems (e.g., 
Bontekoe and van Albada 1987), not representative for a study of binary dynamics, so 
they were avoided in the determination of the final 2D distribution q x L\. The grid 
shown in Table 1 has 24 cells and includes all pairs with q < 6, with 8 cells unoccupied. 

The least massive binary system, (6, 0.5), has a total mass of 2.9 x 10 11 Mq and the 
most massive one, (3,2.0), has a total mass of 1.3 x 1O 12 M0. Such a two-dimensional 
grid defines the mass parameter space of artificial samples generated by Monte Carlo 
simulations. For the simulations of artificial pairs, the 16 occupied cells of Table 1 are 
sampled in proportion to the occupation frequency of the observed sample. 
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Three different orbital eccentricities are considered, namely, e = 0, 0.4 and 0.9, chosen 
as representative of low, medium and high eccentricity orbits. For every simulation panel 
shown in Figure 5, there are 2000 artificial pairs randomly distributed, as described 
above, throughout the 16 {q,L\) bins in Table 1, each pair being generated through the 
following steps, (a) A value of r apo is obtained from the distribution given by eq. 9; (b) 
the value a = a(e) is calculated (eq. 5 or eq. 8, for MOND and MKG, respectively); 
(c) V(r apo ) = a x V r c i rc (r apo ) is calculated; (d) a test is done, using the conservation 
laws of energy and angular momentum, to verify if the pair will ever have a pericentric 
separation (r per ) smaller than 20 kpc (the so-called "merging test", adopted in S90). If 
this is true, the algorithm begins again at (a). Otherwise, it goes to the next step, (e) A 
half-orbit is calculated, i.e., the binary path from apocenter to pericenter. The equations 
of motion, from MOND and MKG potentials, are numerically integrated. Due to the 
characteristics of the force field, this orbit segment is fully representative of the time 
evolution of the system; (f) the half-orbit is rotated by a random angle, in the orbital 
plane. This simple procedure mimics the evolution of the binary over a long time scale. 
(g) An orbital inclination i is obtained according to the distribution F(i) oc sini, i.e., the 
normal to the orbital plane is distributed at random in space, (h) Projected separation 
in the plane of sky (R), and line-of-sight velocity difference (AV) are calculated, at a 
random time instant within the orbit half-period, as well as other relevant quantities, 
such as orbital eccentricity, orbital period, etc. This sequence is the same regardless of 
the gravitational potential adopted. 
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Fig. 5. Monte Carlo simulations of binary galaxies. Left panels represent binaries under MOND 
gravitational potential, and right panels under MKG gravitational potential. The solid curve is 
the fiducial envelope (Keplerian orbits, M/L = 30). The orbital eccentricities are indicated on 
the right upper corner of each panel, and M/L = 5 [in = (L\ + L2) /1Q 11 Lq\. 



Though it is not apparent from Fig. 5, it must be pointed out that, for a given M/L, the 
higher is the total mass of a given pair, the smaller is its AV / L 1 ! 2 ordinate. This does 
not happen with Keplerian models (see eqs. All and A14 ) but both MOND and MKG 
have approximately AV/ L 1 ! 2 oc M -1 / 2 . A similar behavior was also seen with the dark 
halos models worked out in S90. 

Figure 5 shows that MOND and MKG have distinct asymptotic behavior as R in- 
creases. The upper envelope of the MOND circular orbit simulation is independent of R 
while MKG's increases as R 1 ! 2 , as is expected from the linear potential. For non-circular 
orbits, MOND's upper envelopes have an asymptotic value given by a x V c i TC (r apo )/ L 1 / 2 , 
and MKG's again increases with R 1 ! 2 . That is to say, the upper envelope of the MOND 
simulated points will eventually reach a plateau, for large R, and the MKG envelope 
will increase as R 1 ! 2 , since the 1/r component in the potential vanishes in that range of 
separations. The asymptotic limit for the Keplerian counterpart is of course zero. The 
asymptotic behavior of the various models is very promising as a way of discriminat- 
ing amongst them, specially between MOND and MKG. The lack of a sufficiently large 
observed sample of wide pairs is, presently, the main obstacle to such an achievement. 

From Fig. 5, in the range R < 200 kpc, it is readily clear that, for both MOND and 
MKG, binary galaxies in either pure circular orbits or in orbits with only intermediate 
eccentricity (e = 0.4) are not able of explaining the observations since they would require 
extremely high M/L (ps 45, in the range R < 50 kpc) for MOND and MKG standards. 
They also show different trends in the functional dependence of AV / L 1 ! 2 with R, which 
is apparent in the range 50 < R < 100 kpc. Pure high eccentricity orbits (e = 0.9) 
represent viable solutions for both models. Elongated orbits are also consistent with 
the cosmological picture of binary galaxy formation in which such systems are formed 
by galaxies that instead of participating individually of the expanding universe stayed 
together because of their high mutual gravitational attraction, implying returning orbits 
with relatively low angular momentum (e.g., Schweizer 1987). A larger value of a would 
improve MOND's fitting since AV / L 1 ! 2 oc a^ 4 (see eq. 3 for the low acceleration regime). 
In the case of MKG, on the contrary, a smaller value of j would do better since, in the 
large radius range where the linear potential component becomes significant, AV/ L 1 ! 2 
scales pa 7 1 / 2 (see eq. 7), although here a smaller M/L, being still within MKG standards, 
would work as well on lowering the points on the AV / L 1 ! 2 — R plane. It is apparent that 
the only way to distinguish between MOND and MKG lies in the investigation of the 



certainly discriminate between MOND and MKG solutions. 

Although their envelopes are consistent with the observed ones, high eccentricity orbits 
(MOND and MKG panels with e = 0.9) seem inconsistent with the observations because 
the corresponding simulations predict a large concentration of pairs at intermediate and 
large separations, as compared to small ones, and most of observed pairs have small 
separations. Nevertheless, this may be caused by selection effects. It is obvious that a 
combination of binary orbits with eccentricities ranging from small to large values is also 
possible for MOND and MKG. In this case, the high AV/L 1 ! 2 values (with R < 200 kpc) 
will be accounted for by high-eccentricity orbits while the low values will have contribution 
for the whole range of eccentricities. The way to discriminate between the two gravity 
models is an investigation of the relative frequency occupation of the R — AV/L 1 ! 2 plane. 
The large R region (wide pairs), again, is fundamental in this aspect because of the very 
distinct predictions of MOND and MKG in that region (see Fig. 5). To accomplish such 
a discrimination, a full account of sample selection biases and of presence of optical pairs 
must be undertaken. 

The above study gives one a secure indication of the input parameters for a more 
detailed statistical analysis, justifying the initial aim of the work, namely, to present an 
overview of specific applications of alternative theories of gravity to binary galaxy dy- 
namics. The series beginning here will continue with a rigorous analysis of binary galaxy 
samples under the assumptions of MOND models. A third paper is planned on an investi- 
gation of binaries from the point of view of MKG . Both works will extend and complement 
the qualitative study presented here, with a careful account of selection effects on the 
determination of real samples of binary galaxies, with a quantitative measure of sample 
contamination by non-physical pairs, and with the introduction of the same effects in the 
list of prescriptions of Monte Carlo simulations of synthetic samples. 
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The squared relative velocity of two orbiting point-masses Mi and M2 is: 

V 2 (r) = GAl(- - -), (Al) 

where G is the gravitational constant, M = Mi + M2, r is the polar radial coordinate, 
and a is the semi-major axis of the orbit. Consider now a two-dimensional phase space 
represented by the line-of-sight pair relative orbital velocity, normalized by the square 
root of the total pair mass, and the projected pair separation onto the plane of sky. These 
two quantities are related by the projected version of eq. Al: 

/ AV \ 2 

where R is the pair projected separation and AV is the line-of-sight velocity difference, x 
is a projection factor that takes care of the projection of both the orbital relative velocity 
and space separations of the galaxies in the pair: 

AV r G l 1/2 . . 

sin i[cos(f — w) + e cos w], (A3) 



Mi/2 La(l-e 2 ) 
and 

_a(l-e^)_ . n2 . n2 i/2_ 

1 + e cos v ' 

The angles u>, between the line of nodes and the orbit major axis, and i, the orbital plane 
inclination with respect to the plane of sky, give the orientation of the orbit in space. The 
true anomaly (polar angular coordinate) v is the angle, in the orbital plane, between the 
line joining one galaxy to the other and the major axis of the orbit and e is the orbital 
eccentricity. According to eqs. A3 and A4, the projection factor x is given by: 

sin 2 i 2 

X= |cos(i^ — ui) + e cos ui\ x 

1 + e cos v 

1 1 2 

x [1 - sin 2 (i^ - w) sin 2 i] . (A5) 

Note that, for a given eccentricity, x nas a maximum value equal to Xmax(e) = 1 + e 
(the conditions for a maximum projection are u> = 0, i = tt/2 and v = 0. An immediate 
conclusion from eq. A2 is that any observed sample of pure bound pairs would appear 
below the curve 

^G^) 2=xmax(l)G=2G (A6) 

m a AV/M 1 !' 2 - R dia gram. 



that as R approaches zero, AV / M 1 ! 2 becomes increasingly larger. This would be true if 
galaxies were point-masses, which would imply no restriction whatsoever in orbital sizes, 
that is to say, the major axes could have any value. Of course, this is not possible with real 
galaxies; thus, one needs to put a lower limit on orbital sizes. This limit can be defined in 
terms of a minimum allowed pericentric distance, r per = r m . This implies in a minimum 
semi-major axis, a m ; n (e) = r m /(l — e), for a Keplerian orbit with eccentricity e. In fact, 
such condition is equivalent to saying that pairs with pericentric separations smaller 
than r m will merge quickly, which is tantamount to introducing a merging condition on 
Keplerian pairs. 

Upper envelopes, similar to eq. A6, are now derived considering several orbital eccen- 
tricities, and the pericentric cutoff described above. Orbits that have a = a m ; n lead to 
maximum AV / M 1 ! 2 : 

AV r G l 1/2 . . 

sin i[cos(f — w) + e cos w], (A7) 



M 1 / 2 Lr m (l + e) 
at 

D r m(l + e ) f1 -2/ \ • 2 -I 1 / 2 / a o\ 

R=- 1 -— 1 — 8111 \V — LO) 8111 % . (A8) 

(l + ecosi^) L ' 

Substituting the angular conditions for maximum projection (i = 7r/2,w = and v = 0) 
one gets 

AM _lG(l + e)] ^ (Ag) 



that occurs at 

For a given eccentricity there is a well-defined maximum AV/M 1 / 2 ; for all eccentricities, 
maxima occur at the same R = r m . In the case of R > r m the curve that describes the 
maximum values of AV / M 1 ! 2 is 

AV /D , rG(l + e)i 1/2 

and corresponds to orbits with fixed e,a = a m ; n (e) > a m ; n (0) = r m ,i = 7r/2,w = and 
v = 0. In the range < R < r m , one finds the largest possible values of AV/M 1 ! 2 by 
fixing a = a m ; n ,i = 7r/2, and u> = 0, but letting v (that gives the relative position of 
galaxies in orbit) be a free parameter. Doing that in eqs. A7 and A8 one gets: 

AV r G 1 1/2 

(cosi^ + e) (A12) 



M 1 ! 2 Lr m (l + e) 



R: 



»*m(l + e) COS V 



(A13) 



1 + e cos v 

Eliminating v in eqs. A12 and A13 yields the equation of the inner maximum envelope: 

R 



AT/ r G ni/2 

Mi/ 21 ' gJ Lr m (l + e)J Lr m (l + e)-ei? 



+ e 



(A14) 



Note that, for e = (circular orbits), eq. A14 reduces to a line, AV/M 1 ? 2 = (G/r^) 1/2 i?, 
and, for R = r m , i.e., the position of the maximum value of AV/M 1/2 , it takes the form 
of eq. A9, as it should. Figure 6 shows the envelopes for various eccentricities, and for the 
total pair mass given in units of 1O 11 M . The most important point to notice is that the 
merging condition, introduced by only having a > a m ; n , does not allow pairs with high 
AVM 1 ! 2 at small separations (see dashed curves in Fig. 6). 
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Fig. 6. The maximum allowed values of AV/M 1 ^ 2 , as a function of R and eccentricity e, are 
represented by the solid curves (eq. All, for R > r m , and eq. A14, for R < r m ). The dashed 
curves show the envelopes for R < r m , if galaxies are point-masses without restrictions on the 
pericentric separation. The values of the eccentricities are 0, 0.2, 0.4, 0.6 and 0.8. The total pair 
mass is given by Mu = (Mi + M 2 )/W 11 M Q . 
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